Analogy between two-dimensional dislocation systems and layered superconductors 
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We propose a consistent treatment of the applied external stress on a two-dimensional dislocation 
systems by using the analogy between our model and the layered superconductors. Using the results 
of a mathematically rigorous, real-space renormalization-group (RG) calculations in the latter model 
we extend the original model developed by Khantha et al. From our recursion relations a nonlinear 
equation is derived for the transition temperature as a function of the applied stress. Our results are 
compared with the predictions of Khantha's model for three particular materials and the reasons 
for the deviations are discussed. Possible extensions of our model are outlined. 
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In the context of dislocation pattern formation the col- 
lective behavior of the dislocation systems seems to play 
a fundamental role and a new approach has been used 
to determine the distribution function of the dislocation 
systems [Q. However, in this model the temperature de- 
pendence has not been included. Recently, a new type 
of mechanism of the cooperative generation of the dislo- 
cations has been proposed by Khantha, Pope and Vitek 
(hereafter referred as to KPV model) §- j§. The KPV 
model of the unstable generation of dislocations driven 
by thermal fluctuation and aided by an applied stress 
on the solid is based on the theory of the screening of 
topological defects (vortices) developed by Kosterlitz and 
Thouless [QH) and Berezinskii M (hereafter referred as to 
KTB). The KTB type transition of the topological de- 
fects (vortices) has already been well established in sev- 
eral physical systems such as dislocation mediated melt- 
ing in 2D crystal [ l(],|ll| , layered high-temperature super- 
conductors [p^p^ l, Coulomb gas mM, superfluid 4 Helium 
films, 2D X-Y model, liquid crystals (for review see {Hi!). 
The KTB phase transition seems to play a fundamental 
role in the behavior of some mesoscopic devices like dou- 
ble layer quantum Hall systems Jl6| , where the statistical 
mechanics of the so called merons analogous to that of 
the vortices. 

In the KTB type transition the bounded dislocation 
pairs unbind (free dislocations appear in thermal equilib- 
rium) at a certain temperature where an universal jump 
in shear modulus takes place. As an extension of the 
KTB theory in the KPV model the effect of applied stress 
is included additionally through a mean-field type calcu- 
lation. In the model one of the main results is the recur- 
sion equations for the coupling constant and the fugacity 
which relate to the strength of interaction between dislo- 
cations and the average number of dislocations in the sys- 
tem, respectively. Due to the applied external stress the 
change in fugacity was taken into account in the deriva- 
tion of the recursion relations and it results in a modified 
recursion equation for the fugacity in comparison with 
that of the KTB theory. 

However, the stress not only increases the fugacity (and 



consequently the density of dislocation pairs) but indi- 
rectly affects the coupling constant which can be un- 
derstood as follows. The number of dislocation pairs 
increases with increasing stress. Larger number of dis- 
location pairs makes more effective the screening of the 
interaction between dislocations resulting in a decrease 
in the coupling constant. It is clear that in a consistent 
treatment of the effect of the stress one needs to take 
into account this indirect effect of the stress which was 
omitted in the KPV model. 

The vortex gas in the layered superconductors shows 
many common aspects with the above discussed disloca- 
tion systems in the sense that the character of the inter- 
action between vortices and dislocations can be identical 
in some cases. The Hamiltonian used in the KPV model 
is the same as that in the study of the high temperature 
layered superconductors by taking the 2D limit of the 
work of [jl3| for the case of a vortex gas with a current 
applied uniformly through all layers. This current can 
be associated with the applied stress in the dislocation 
systems. Comparing the Hamiltonian of the two models 
one can identify the corresponding parameters between 
them. This identification makes possible to obtain the 
recursion relations for dislocation systems by using the 
recursion relations found in the case of layered supercon- 
ductors in the 2D limit (see Pierson's paper in Ref . |0 ) . 
Pierson in his work performed a mathematically rigorous, 
real-space renormalization group (RG) analysis following 
the original approach used first by Kosterlitz || . 

Of course, one can derive the recursion equations inde- 
pendently of using the above analogy. Since in Pierson's 
work the derivation of the recursion relations are not pre- 
sented (only the results) we performed the RG analysis 
with our dislocation systems for clarity. We followed the 
procedure used in Ref. [ll] (the details to be published 
somewhere else). 

In this letter we present a consistent treatment of the 
dislocation systems in the presence of applied stress. We 
use the same model as in Ref. g (Eq. (1) therein) and 
give the new recursion relations which includes the effect 
of the applied external stress both in the fugacity (as in 
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the KPV model) and in the coupling constant. In our 
calculations the stress is taken in lowest nonvanishing or- 
der both in the equation of the fugacity and the coupling 
constant. From this study it turns out that the equation 
for the fugacity remains unchanged (same as the second 
equation of (4) in Ref. Q), while in the equation for the 
coupling constant an additional term proportional to the 
square of the stress appears. Therefore, our calculations 
reveal that the indirect effect of the stress modifies only 
the recursion relation for the coupling constant. How- 
ever, this latter change has several consequences which 
will be discussed below. 

In the KPV model the stress dependence of the tran- 
sition temperature is given by a nonlinear implicit equa- 
tion derived from the stability analysis of the recursion 
equations (see Eq. (6) in Ref. ||). Using our new re- 
cursion relations a different nonlinear equation is found 
for the stress dependence of the transition temperature. 
The above mentioned indirect effect of the stress can be 
associated precisely. Then, for different materials, we 
compare the results obtained from the KPV model with 
that found by using our recursion relations. Qualitative 
explanations are given for the deviations of the transition 
temperature from the KPV model. Finally, we summa- 
rize the necessary steps for improving the model in order 
to treat dislocation systems more realistically. 

A consistent treatment of the applied stress and the 
role of the indirect effect of the stress in the dislocation 
systems are our main subject in this paper. Our results 
obtained from the extension of the KPV model are a new 
contribution in this field. 

In the frame work of the continuum elasticity theory 
p7| using the same interaction energy between disloca- 
tions as in the KPV model (Eq. (Al) in Ref. g) the 
Hamiltonian of the system Ti. containing N number of 
dislocations may be expressed as 

N 
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where j3 — l/k^T and ry = — rj is the distance be- 
tween the ith and jth dislocation. The summation over 
pairs < i,j > in the Hamiltonian assumes i ^ j and 
counts each pair just once. The charge pi = +1 if the 
Burgers vector at site i is directed to the positive x-axis 
and equals to —1 if it is pointed to the negative x-axis. 
The charges satisfy the neutrality condition pi = 0. 

The interactions between dislocations is given by the 
first term in ([!]) where 

G(r)=ln(|r|/oo) (2) 
and clq is the cut-off radius. The coupling constant in (jl|) 
K Q = H and J = ££+g where ^ and B arc the 



shear and bulk moduli in the absence of dislocations and 
bo is the magnitude of the Burgers vector. 

The second term in the Hamiltonian is a sum of the 
interactions between dislocations and the applied exter- 
nal stress that can be derived from the Peach-Kochlcr 
formula This interaction may be written in such 

a form that is equivalent to the interaction energy p^Er; 
between the 'charge' pi and the homogeneous external 
'electric field' E = a^box, where x is the unit vector 
pointed to the positive x-axis. Then Eq in the Hamilto- 
nian is related to the applied external stress by 

E =a \E\ =a/3b a (3) 

and the interaction potential W is given by 

W/(r 2 ) = -r !X * (4) 
a 2 

Note that a constant —1/2 has been added to W to be 
consistent with the definition of the core energy E c of 
the dislocation (which must be one half the energy of a 
dislocation pair at smallest separation and oriented along 
the direction of the external field E) . 

Finally, in the last term of the Hamiltonian, the fugac- 
ity yo is related to the core energy by 

\ny = -/3E c + E /2 7 (5) 

where the second term is due to the shift of the potential 
W. 

It is easy to show that our Hamiltonian is equivalent to 
that of the KPV model and in this form one can also rec- 
ognize the resemblance to the Coulomb gas model with 
applied external electric field. Furthermore, it is clear 
that the Hamiltonian ([j]) is similar to that of the model 
of the layered superconductors if we omit the terms con- 
taining A which relates to the interlayer coupling (see Eq. 
(1) - (3) in Ref. fl3|l ). Below we make a precise corre- 
spondence between them. 

First, as an extension of the KPV model we derive the 
recursion equations including the effect of the stress both 
in the equation of the fugacity and the coupling constant. 

A possible starting point for such a calculation is the 
partition function. Following the procedure used in Ref. 
|TTf one can derive the recursion relations (the details to 
be published somewhere else). 

However, using the above discussed analogy between 
our dislocation systems and the layered superconductors 
we may obtain our recursion relations by simple iden- 
tifying the parameters in the two models via the two 
Hamiltonians. We found the following correspondences 
between the two models (the parameters in Ref. [jD| are 
on the left hand side, while our parameters are on the 
right hand side) 

r^a, e -*• I, (3p 2 -> 2nK, (6) 
2 

x->—-l, y^2ny, J (3rpi -> Ep i: (7) 
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where J does not contain the absorbed factor of \[0t as 
in Ref. pl| after Eq. (7). Note that p 2 — 1 in our case. 

Using the above identifications and the recursion rela- 
tions given in Ref. fl3|] and omitting the A terms we find 
for the coupling constant K, the fugacity y and E 
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Eqs. (P)-(|To|) are our final coupled differential equa- 
tions with the initial conditions K(l = 0) = K , E(l = 
0) = E Q and y(l = 0) = enp{-f3E c + E Q /2). 

The indirect effect of the stress in the coupling constant 
is reflected by the second term in (^). Omitting the E 2 /4 
term in Eq. (||) we recover the same recursion equations 
as in the KPV model. The equation for the fugacity given 
in (^|) is the same as in the KPV model. Therefore, the 
inclusion of the indirect effect of the stress changes only 
the renormalization of the coupling constant K. It is 
also clear that this indirect effect in K appears in second 
order in a since E ~ a. 

At fixed stress the flows obtained from our recursion 
relations show that below at certain temperature T c the 
fugacity y — > oo and for T > T c y — ► 0. The fugacity 
y is related to the mean density of the unbound dislo- 
cation pairs in the system. Increasing the applied stress 
(note that Eq ~ a) T c decreases. According to our re- 
cursion relations the reason is threefold: (i) the stress 
increases the fugacity yo, (ii) the last term in Eq. ([)]) has 
positive feedback for y, (iii) in Eq. (||) the coupling con- 
stant K lowers, i.e., the dislocation pairs are more weakly 
bounded. This latter effect, which is our central issue in 
this letter, is related to the above mentioned indirect ef- 
fect of the applied stress. 

In the KPV model a nonlinear equation was derived 
for the stress dependence of the transition tempera- 
ture by linearizing the recursion equations around the 
fixed point. Similar way, using Eqs. (|^)-([Io|) we have 
found that the dimensionless transition temperature t c = 
k-QTc/ J is a solution of the following nonlinear equation 



t, = 



1 - 2ttC(E ) exp (-0E C + E /2) 



En 



(11) 



where C(E ) 
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o-f is the dimensionless stress while V is the 



and a' 

magnitude of the Burgers vector in units of an. In the 
KPV model C(E ) = 1. The indirect effect of the applied 
stress in Eq. ( pi] ) appears via the E 2 term in C(En). The 
expressions for t c obtained from the two models differ 
only by the parameter C{Eq). 



To study the role of the indirect effect of the stress as 
an extension of the KPV model, we solved numerically 
Eq. (11) for the transition temperature t c at different 



values of a. We have chosen three particular materials, 
namely, TiAl, Si and NiAl with the same material pa- 
rameters as in 0]. According to these parameters the 
core energies E c = 0.25J for TiAl, E c = 0.78J for Si 
and E c = 0.27 J for NiAl were used. Fig. [l] shows the 
variation of t c with the applied stress, a obtained from 
our calculations and from the KPV model for TiAl and 
Si. For comparison, the transition temperature is plotted 
in unit of J/fce, while we used the same stress range for 
both materials. 
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FIG. 1. The transition temperature, t c as a function of the 
applied stress, a for TiAl and Si. For parameters see the text. 
tc and the stress are units of J/fce and GPa, respectively. In 
case of TiAl the solid line and the line with circles correspond 
to the KPV and our model, respectively. In case of Si the 
dashed line and the line with squares are the results from the 
KPV and our model, respectively. 

The result from the KPV model can be calculated by 
simply taking C(Eq) — 1 in Eq. (pi]). It is important 
to mention that in both the KPV and our model it is 
assumed that the fugacity, i/ 2 C 1, i.e., the average den- 
sity of dislocations is small enough. From our calcula- 
tion it turns out that for a < 3.5 GPa y 2 < 1. For all 
the three materials the fugacity was always lower in our 
model comparing with the KPV model. 

From Fig. [l] one can see that the indirect effect of the 
stress results in a small deviation from the original KPV 
model for Si, especially for small values of stress. In- 
creasing the stress the differences become larger but it 
is still about 7% at a = 3 GPa. However, the situa- 
tion for TiAl is changed substantially. Qualitatively, two 
changes can be seen from the figure. First, the transition 
temperature is smaller for TiAl than for Si for all values 
of stress (it is approximately 0.1 smaller at zero stress). 
Second, the transition temperature calculated from our 
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model departs from the prediction of the KPV model at 
lower values of a for TiAl comparing with those for Si. 
The differences in t c between our and the KPV model 
are substantially larger for TiAl than in the case of Si 
(the difference is about 25 % at u = 3 GPa for TiAl, 
while only 7% for Si). We have also calculated the stress 
dependence of t c for NiAl. We found that the qualitative 
feature remains the same as for TiAl and only a slight 
change can be seen in the values of t c vs. a (for zero stress 
t c = 0.112 for NiAl). 

The first observation can be understood from the core 



energy dependence of t c at zero stress, using Eq. (11) 



With increasing E c the transition temperature decreases 
due to the exponential term in Eq. ( |ll| ) . The physical ex- 
planation of this behavior is as follows. For large E c the 
number of thermally activated dislocations is relatively 
small. This implies that the screening of the disloca- 
tion pairs by other dislocations is weaker. The coupling- 
strength (normalized by J) between pairs remains strong. 
Hence, the unbounding of the pairs, which is related to 
the phase transition, takes place at higher temperature. 
The core energies for TiAl and NiAl are approximately 
the same, while for Si it is much larger. Thus, t c will be 
higher for Si than TiAl and for NiAl it is nearly the same 
as that for TiAl. 

The second observation is less obvious. Again, the 
analysis of the core energy dependence of t c shows that 
at lower E c the indirect effect of the stress becomes 
more prominent. The parameter C(Eq) in Eq. ( pT| ) can 
be dominant besides the suppression of the exponential 
term. On the other hand, with inrcasing stress, t c ob- 
tained from our model is always lower than that predicted 
from the KPV model independently from the value of 
E c . The physical reason is that at higher a the so-called 
'electric field' makes weaker the bounds between the dis- 
location pairs. In Eq. (ph it is reflected by the second 
term. This term is entirely the consequence of the indi- 
rect effect of the stress. 

In the KPV model the Burgers vector has only two 
possible orientations correspondig to +1 and —1 charges. 
However, in a solid there can be more orientations of the 
Burgers vector depending on the crystal structure. It is 
desirable to extend the model in order to take into ac- 
count the possible Burgers vector systems in the crystal. 
In the present study (and also in the KPV model) the 
role of the angular dependent force between dislocations 
has not been considered. The quantative behavior of the 
stress dependence of the transition temperature is likely 
to influence by this force. Some preliminary work on 
these directions have already been made pQ| . 

The KPV model has been applied to the brittle-to- 
ductilc transition and the deformation of whiskers in Ref. 
||. We believe that our model contributes to a better 
understanding of these two phenomena. 

In conclusion, we present an extension of the KPV 
model by taking into account the indirect effect of the 



applied external stress on the solid. Realizing the anal- 
ogy between our dislocation systems and the layered su- 
perconductors we derive the correspondig recursion equa- 
tions. From these equations we obtain a nonlinear equa- 
tions for the transition temperature. The implication of 
our model is discussed in the case of three particular ma- 
terials. In generaly, we find that the indirect effect of the 
stress in T c plays crucial role in some materials depending 
on their core energies. 
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